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Abstract 

We present the results of a time-dependent Dalitz plot analysis of — 7r+7r~7r° decays based on 
a 414fb^^ data sample that contains 449 X IQ^BB pairs collected on the T(4S') resonance with the 
Belle detector at the KEKB asymmetric energy e'^e~ collider. The direct CP violating parameters 
of B^ p^-K^ decay mode are measured to be = +0.22±0.08(stat.)±0.05(syst.) and = 
+0.08±0.17(stat.) ib0.12(syst.). We also measure the CP violating parameters of the decay mode 
_^ p^ttO asApO^o = -0.45±0.35(stat.)±0.32(syst.) and<Sp0^o = +0.15±0.57(stat.)±0.43(syst.). 
Combining our analysis with information on charged B decay modes, we perform a full Dalitz and 
isospin analysis and obtain a constraint on the CKM angle 02, 4>2 = (83l23)°- ^ large CKM- 
disfavored region (^2 < 8° and 129° < (P2) also remains at 68.3% confidence level. 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 
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1. INTRODUCTION 



In the standard model (SM), CP violation arises from an irreducible phase in the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix [1]. The SM predicts that measurement of a 
CP asymmetry in time-dependent decay rates of and gives access to the CP violating 
phase in the CKM matrix [2-4]. The angle 02 of the CKM unitarity triangle can be mea- 
sured via the tree diagram contribution in 6 — >^ uud decay processes, such as 5° — > 7r"'"7r~, 
B^ — > p^TT^, or B^ — > p^p~ [5]. In these decay processes, however, contributions from 
so-called penguin diagrams could contaminate the measurement of 02- Snyder and Quinn 
pointed out that a Dalitz plot analysis of B^ pn, which includes B^ P'^'^~ , B^ P~'^~^, 
and 5° p^TT*^, offers a unique way to determine 02 without ambiguity. The Dalitz plot 
analysis takes into account a possible contamination from a 6 — > (i penguin contribution [6] . 
In addition, an isospin analysis [7, 8] involving the charged decay modes, S+ p+7r° and 
S+ — > p^TT"*", provides further improvement of the 02 determination. 

1-1. KEKB and Belle Detector 

KEKB [9] operates at the T(45') resonance {^/s = 10.58 GeV) with a peak luminosity 
that exceeds 1.6 x 10^^ cm^^s"^. 

At KEKB, the T(45') is produced with a Lorentz boost of — 0.425 nearly along the 
electron beamlinc (z). Since the 5° and B^ mesons are approximately at rest in the T(4jS') 
center-of-mass system (cms), At can be determined from the displacement in the z direction, 
Az, between the vertices of the two B mesons: At ~ Az/(3jc. 

The Belle detector is a large-sohd-angle magnetic spectrometer that consists of a sihcon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cercnkov counters (ACC), a barrel-like arrangement of time-of-fiight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of Csl(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect Kl mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [10]. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
data sample of 152 x 10^ BB pairs (DS-I), while a 1.5 cm beampipe, a 4-layer silicon detector 
and a small-cell inner drift chamber were used to record the remaining 297 x lO^BB pairs 
(DS-II) [11]. 

1-2. Analysis flow 

The analysis proceeds in the following steps. First, we extract the signal fraction (Sec. 2). 
We then determine the sizes and phases of the contributions from radial excitations (Sec. 3). 
Using the parameters determined in the steps above, we perform a time-dependent Dalitz 
plot analysis (Sec. 4-6). The fit results are interpreted as quasi-two-body CP violation 
parameters (Sec. 7) and are used to constrain the CKM angle 02 (Sec. 8). 
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1-3. Differential decay width of time-dependent Dalitz plot 



We measure the decay process 5° — > n'^n'n^, where we denote the four-momenta of 
the TT"*", 7T~ , and tt^ by p+, p_, and po, respectively. The invariant- mass squared of their 
combinations 

s+ = (p++po)^ s_ = (p_+po)^ so^(p++p-f (1) 

satisfies the following equation 

s+ + s_ + So ^ m|o + 2m^+ + m^o (2) 

by energy and momentum conservation. The differential (time-integrated) decay width with 
respect to the variables above (Dalitz plot) is 

where A^^ is the Lorentz-invariant amplitude of the S°(S°) — > 7r+7r~7r° decay. 

In the decay chain T{AS) B'^B^ fif2, where one of the -B's decays into final state 
fi at time ti and the other decays into another final state /2 at time t2, the time dependent 
amplitude is 

^(^1,^2) ~ e-(r/^+^^)(*^+*^){cos[Am,(ii-i2)/2](>liA2-Xi^) 

-z sm[Amd{h - t2)/2] (^^1^2 - ^ ^^2) } ■ (4) 

Here, p and q define the mass eigenstates of neutral B mesons as pB^ ± with average 
mass M and width F, and mass difference Am^. The width difference is assumed to be zero. 
The decay amphtudes are defined as follows, 

(5) 
(6) 
(7) 
(8) 

In this analysis, wc take A^t^ as Ai and choose /2 to be a flavor eigenstate, i.e., = or 
A2 = 0. Here we call the B decaying into /i = fcp = n~^7i~n^ the CP side B while the 
other B is the tag side B, f2 — /tag(/tag)- The differential decay width dependence on time 
difference At = tcp — ^tag is then. 



Ai 


= A{B^ - 




Ai 


= A(W - 




A2 


= A{B° - 


-/2 


M 


= A(W - 


^f2 



dr^e-ri^*i{(|A3.|2 + |A3 

-gtag • (|^3^P - |^3^n cos(AmdAi) + gtag • 2Im [jA*^^3^J sin (Am^Ai) jo? Ai , 

(9) 

where we assume \q/p\ = 1 {CP and CPT conservation in the mixing) and |/l(i?° — >■ /tag)| = 
|A(i?° —>■ /tag)L ^iid integrate over tsum = tcp + ttag- Here gftag is the 6-flavor charge and 
?tag — +1(~1) when the tag side B decays into a B^ (B^) flavor eigenstate. 
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Combining the Dalitz plot decay width (3) and the time dependent decay width (9), we 
obtain the time dependent Dahtz plot decay width 



dr ~ \A{At; s+, dAt ds+ ds_ , (10) 

where 



P) cos(Am<iAi) + 5tag • 2Im [m*^A3^] sin(AmdAi)} , 



(11) 



^37r = A3^{S+, S_) , As^ = ^37r(s+, S_) . (12) 

Although there could exist contributions from decays into non-p7r tt^tt'tt^ final states, 
such as /o(980)7r°, /o(600)7r'^, cj7r°, and non-resonant 7r"'"7r~7r°, we confirm that these contri- 
butions are small and their efi^ects are taken into account as systematic uncertainties (Sec. 6 6- 
5). We therefore assume that the — > 7r+7r~7r° decay is dominated by the S° — (p7r)° 
amplitudes: S° P+tt", S° p~7r+, and S° p°7r°, where p can be p(770), p(1450), or 
p(1700). The Dalitz plot amphtude ^43^(5+, s_) can then be written as 



^3.(s+, s.) = /+(s+, s.)A+ + /_(s+, s.)A- + /o(s+, s_)AO , (13) 
^As^is+,s.) = J+is+,s.)A^ + J_{s+,s.)A- + Us+,s.)A' , (14) 

where the functions /k(s_|_, s_) (with p charge k — +,—, 0) incorporate the kinematic and 
dynamical properties of the B^ decay into a vector p and a pseudoscalar tt. The goal of this 
analysis is to measure the complex coefficients ^4+, A~ , A^, A'^, A~, and A^; and constrain 
the CKM angle 02 using them. 



1-4. Kinematics of — > (P^)'' 

The function /k(s+, s_) can be factorized into two parts as 

y,(s+, s_) = T^F^is,) {k = +, -, 0) , (15) 

where FTr{sK) and Tj correspond to the lineshape of p and the helicity distribution of the p, 
respectively. Here we assume that a single unique functional form for the lineshape F^(s) 

can be used for all six '/^[27]. Since this assumption has no good theoretical or experimental 
foundation, we check the goodness of the assumption with data and assign systematic errors. 

The lineshape is parameterized with Breit-Wigner functions corresponding to the p(770), 
p(1450), and p(1700) resonances: 

F^{s) = BWp^rro) + P ■ BW,^U50) + 7 • 5W^p(i7oo) , (16) 

where the amplitudes /3 and 7 (denoting the relative size of two resonances) are complex 
numbers. We use the Gounaris-Sakurai(GS) model [12] for the Breit-Wigner shape. 
In the case of pseudoscalar- vector {J — 1) decay, Tj is given by 

7^r = -4|p,||Pfe|cos^^'\ (17) 
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T+ = -4|pV||p_| cos^+- , \ 

Tf = -4|po I |p+ 1 cos ^^0+ , (18) 
ro = -4|p_||po|cosro, / 

where Pj-,Pk are the three momenta of the vr^ and tt'^ in the rest frame of p'^ (or the ttV'' 
system), and the 0^^[= 0^) is the angle between pj and pk (see Fig. 1). 



i 




FIG. 1: The relation between three pions in the rest frame of p'^. 



1-5. Fitting parameters 

Inserting (15), into expressions (13) and (14), the coefficients of equation (11) become 
Ke{+,-,o} /t<CTe{+,-,o} 



^ ^ Ke{ + ,-,0} K«TG{ + ,-,0} 



(20) 



with 



= (lA'^p ± |A''|2)/iv , (21) 
^±,Re(im) _ Re(Im) [AM"* ± M"*] /iV , (22) 

4 = Im [A'' A''*] /N , (23) 
J^; = Re p'^A'^* - 'A'A''*] /N , (24) 
/^^ = Im p'^A^* + A^A""*] /N , (25) 

where A'^ is a normahzation factor. The 27 coefficients (21)-(25) are the parameters deter- 
mined by the fit. This parameterization allows us to describe the differential decay width 

as a linear combination of independent functions, whose coefficients arc the fit parameters 
in a well behaved fit. We fix the overall normalization by requiring = 1, i.e., we use 
= \A~^\'^ + I^^P- Thus, 26 of the 27 coefficients are free parameters in the fit. 
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1-6. Square Dalitz plot (SDP) 



The signal and the continuum background e+e" — > qq{q = u,d,s,c), which is the dom- 
inant background in this analysis, populate the kinematic boundaries of the usual Dalitz 
plot as shown in Figs. 2 and 3. Since we model part of the Dalitz plot probability density 
function (PDF) with a binned histogram, the distribution concentrated in a narrow region 
is not easy to treat. We therefore apply the transformation 

ds+ds- — > I det J\dm'd9' , (26) 

which defines the square Dalitz plot (SDP) [13]. The new coordinates are 

m'= - arccos 2 — L - 1 , (27 

\ ™max ™min ; ' v / 



TT 



"0 



0' ^ -9o -e-A , (28) 

TT \ TT / 

where mo — ^/so, m^^^ — m^o — 171^^0 and m™"^ = 2m^+ are the kinematic limits of mo, and 
J is the Jacobian of the transformation. The determinant of the Jacobian is given by 

I det J| = 4|p+||po|mo • — —nsm{7rm') ■nsm{ne') , (29) 

where p+ and po are the three momenta of 7r+ and 7r° in the tt+tt" rest frame. 
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FIG. 2: Distribution of signal Monte Carlo (without detector efHciency and smearing) in the Dalitz 
plot. 
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(b) Square Dalitz plot 



FIG. 3: Distribution of qq background (from the data Mbc sideband) in the Dalitz plot. 



9 



2. EVENT SELECTION AND RECONSTRUCTION 



To reconstruct candidate 5° tt+tt^tt^ decays, charged tracks reconstructed with the 
CDC and SVD are required to originate from the interaction point (IP) and to have trans- 
verse momenta greater than 0.1 GeV/c. We distinguish charged kaons from pions based on 
a kaon (pion) hkehhood Ck(-w) derived from the TOF, ACC and dE/dx measurements in 
the CDC. Tracks that are positively identified as electrons are rejected. 

Photons are identified as isolated ECL clusters that are not matched to any charged 
track. We reconstruct 7r° candidates from pairs of photons detected in the barrel (end-cap) 
ECL with > 0.05 (0.1) GeV, where E^ is the photon energy measured with the ECL. 
Photon pairs with momenta greater than 0.1 GcV/c in the laboratory frame and with an 
invariant mass between 0.1178 GeV/c^ to 0.1502 GeV/c^, roughly corresponding to ±3cr of 
the mass resolution, are used as 7r° candidates. 

We identify B meson decays using the energy difference AE = E^^ — E'b™^ 
beam-energy constrained mass Mbc = ^/(-E'b^m)^ ^ (Pb^*^)^) where i?bcam beam energy 

in the cms, and E'^'^^ and p^** are the cms energy and momentum of the reconstructed B 
candidate, respectively. 

We select candidates in a fit region that is defined as —0.2 GeV < AE < 0.2 GeV and 
5.2GeV/c^ < ilfbc < 5.3GeV/c^. The fit region consists of the signal region that is defined 
as -0.1 GeV < AE < 0.08 GeV and Mbc > 5.27GeV/c^ and the complement, called the 
sideband region, which is dominated by background events. 

The vertex position for the B^ — > 7r+7r~7r° decay is reconstructed using charged tracks 
that have enough SVD hits [14]. The /tag vertex is obtained with well reconstructed tracks 
that are not assigned to fcp- A constraint on the interaction-region profile in the plane 
perpendicular to the beam axis is also used with the selected tracks. 

The 6-fiavor of the accompanying B meson is identified from inclusive properties of parti- 
cles that are not associated with the reconstructed B^ fcp decay. We use two parameters, 
the 6-flavor charge gtag and r, to represent the tagging information [15]. The parameter r is 
an event-by-event, Monte Carlo (MC) determined flavor-tagging dilution factor that ranges 
from r = for no flavor discrimination to r = 1 for unambiguous flavor assignment. It 
is used only to sort data into six r intervals. The wrong tag fractions for the six r inter- 
vals, wi (/ = 1,6), and differences between 5° and B^ decays, Awi, are determined using a 
high-statistics control sample of semileptonic and hadronic b ^ c decays [15-17]. 

The dominant background for the B^ — > 7r"'"7r~7r'' signal comes from continuum events. 
To distinguish these jet-hke events from the spherical B decay signal events, we combine a 
set of variables that characterize the event topology into a signal (background) hkelihood 
variable £sig(bkg), and impose requirements on the likelihood ratio T^-s/b = >Csig/(>Csig -|- >Cbkg) 
that depend on the quality of flavor tagging. 

When more than one candidate in the same event is found in the flt region, we select the 
best candidate based on the reconstructed 7r° mass and T^-g/b- 

After the best candidate selection, we reconstruct the Dalitz variables s+, sq and s_ from 
1) the four momenta of the vr"*" and vr", 2) the hclicity angle of the (i.e., the helicity 
angle of the tt'^ti" system), and 3) the relation: m^o^ + 2m^+^ + m^o^ = s+ -|- s_ -|- Sq- 
Note that the energy of 7r° is not explicitly used here, which improves the resolution of the 
Dahtz plot variables. We reject candidates that are located in one of the following regions 
in the Dalitz plot: ^ > 0.95 GeV/c^ and ^ > 1.0 GeV/c^ and ^ > 1.0 GeV/c^; 
^/s^ < 0.55 GeV/c^ or < 0.55 GeV/c^ or < 0.55 GeV/c^. In these regions, the 
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fraction of pn signal is small. However, radial excitations (p(1450) and p(1700)) are 

the dominant 7r"'"7r~'7r° contributions in the region with > 1.0 GeV/c^. Since 

the amplitudes of the radial excitations are in general independent of the amplitude of the 
p(770), they are considered to be background in our analysis; vetoing the high mass region 
considerably reduces the systematic uncertainties due to their contributions. 

Figure 4 shows the Mbc and AE distributions for the reconstructed B^ —>■ tt'^tt'tc^ can- 
didates within the AE and Mbc signal regions, respectively. The signal yield is deter- 
mined from an unbinned four- dimensional extended-maximum-likelihood fit to the AE- 
Mbc and Dalitz plot distribution in the fit region defined as Mbc > 5.2GeV/c^ and 
—0.2 GeV < AE < 0.2 GeV, where the Dalitz plot distribution is only used for the events 
inside the AE-M]^^ signal region. The AE-M]^^ distribution of signal is modeled with binned 
histograms obtained from MC, where the correlation between AE and Mbc, the dependence 
on PttO; and the difference between data and MC are taken into account. We also take 
into account incorrectly reconstructed signal events, which we call self cross feed (SCF) and 
which amounts to ~ 20% of the signal. In a SCF event, one of the three pions in fcp is 
swapped with a pion in /tag or the 7r° in fcp is misreconstructed. We give the details of 
the Ai?-Mbc and Dalitz plot PDFs of the SCF events in appendix A. For continuum, we 
use the ARGUS parameterization [18] for Mbc and a linear function for AE. The AE-M-^c 
distribution of BB background is modeled by binned histograms based on MC. The Dalitz 
plot distributions for all components are modeled in the same way as the time- dependent fit 
described later, but integrated over the dimensions of the proper time difference. At, and 
the flavor of the tag side B, gtag- The fit yields 987 ± 42 5° ^ tt+tt^tt^ events in the signal 
region, where the errors are statistical only. 
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FIG. 4: The Mbc (upper) and A£'(lower) distributions within the AE and Mbc signal regions. The 
histograms are cumulative. Solid, dot-dashed, dotted and dashed hatched histograms correspond 
to correctly reconstructed signal, SCF, BB, and continuum PDFs, respectively. 
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3. DETERMINATION OF THE CONTRIBUTIONS FROM RADIAL EXCITA- 
TIONS 



Using the same data sample as described above but performing a time-integrated 
Dalitz plot fit with a wider Dalitz plot acceptance, 0.55GeV/c^ < ^/sq < 1.5GeV/c^ or 
y/s^ < 1.5 GeV/c^ or ^^sZ < 1.5 GeV/c^, we determine the p lineshape, i.e. the phases and 
amplitudes of the coefficients (3 and 7 in equation (16), which we use for all of the decay 
amplitudes. In this fit, we use the PDG values [19] for the masses and widths of the p(1450) 
and p(1700). The fit yields 

1/31 = 0.30^°:°^ arg/5=(213l}^r, I7I = 0.07 ± 0.03 , arg7 = (911^)° . (30) 

The mass distributions and fit results are shown in Fig. 5. Figure 6 schematically shows 
how the radial excitations contribute to our fit result. Note that the values given above for 
(3 and 7 and their errors are not meaningful measurements of physics parameters but rather 
are quantities needed for the time-dependent Dalitz fit. This is because these parameters 
are determined from the interference region, the interference between p+vr^ and p^ir^, etc., 
and depend on the unfounded common lineshape assumption of the equation (15). However, 
because statistics are still low, the time-dependent Dalitz analysis would not be possible if 
we were to discard the common lineshape assumption. 

Thus, it is important to determine the common or average lineshape as well as obtain an 
upper limit on the deviation from the average lineshape for each of the 6 decay amplitudes. 
For this purpose, we put constraints on additional amplitudes that describe 1) the excess in 
the high mass region, ^/s > 0.9GeV/c^, and 2) the interferences between radial excitations 
and the lowest resonance, the p(770): interferences between p(770)"'"7r~ and p(1450)~7r"'", 
etc. The nominal fit is performed with the average lineshape determined above, fixing 
all of the additional amplitudes to zero. When fioating the additional amplitudes for the 
other resonances, we obtain results consistent with zero for all of the additional amplitudes 
but with large uncertainties compared to the errors for the average lineshape parameters 
above. We use the fit result with the additional lineshape parameters fioating including their 
uncertainties in the systematic error study. 




FIG. 5: Mass distributions and fitted lineshapes in p'^ir" (left), p^vr^ (middle), and p^ir^ (right) 
enhanced regions. The histograms are cumulative. Solid, dot-dashed, dotted and dashed hatched 
histograms correspond to correctly reconstructed signal, SCF, BB, and continuum PDFs, respec- 
tively. Note that there are feed-downs from other quasi-two-body components than those of interest, 
especially in the high-mass regions. For example, the high-mass region (mo ^ l.OGeV/c^) of the 
enhanced region (right) includes large contributions from p^yr^. 
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FIG. 6: A schematic figure of the fit result of the hneshape and the contributions from radial 
excitations. Note that our definition of i^7r(s) does not include the factor 1/(1 + /? + 7) as in 
equation (16). One can see that the /o(770) and /9(1450) destructively interfere with each other 
near ^/s = itIt^t^ = 1.4(GeV/c^), which means that the p{1450) has a large impact on the phase of 
F7r(s) although the absolute value of |-F7r(s)| is not much affected. 



4. TIME-DEPENDENT DALITZ PLOT ANALYSIS 

To determine the 26 Dalitz plot parameters, we define the following event-by-event PDF: 

Pix) = /sig j9sig(AE,Mbc;m',6''; A)f:,gtag;p,ro) 
+fBB PBsi^E, Mbc; m', 6'- At, gtag) 

p[g{AE, Mbe; m', 9'; At, gtag) , (31) 

where Pgjg, p''^^ and p^^ are PDF's for signal, BB background and continuum background 
in flavor tagging region /, respectively, and /sig, fsB and fqq are the corresponding fractions 
that satisfy 

/sig + IbB + /<jg = 1 • (32) 
The vector x, the argument of p, corresponds to a set of event-by-event variables: 

X = {AE, Mbe; m', 9'; At, gtag, I; P.o) . (33) 

The PDF for signal events consists of a PDF for the correctly reconstructed events p[j.^c 
PDFs for SCF(Self Cross Feed) events p[ {i = NR, CR): 

^ pI^,{AE, Mbc; m', 9'; At, gtag;p.o) + J2i=NRCRPii^E^ ^bc; m', 9'; At, gtag) 
Psig = 7 — , (34) 

^truc + 2^j=NR,CR 

where NR and CR represent the 7r° (neutral) replaced and (charged) replaced SCF, 
respectively, and n are the integrals of the PDFs. A detailed description of the PDF for 
each component can be found in appendix A. 
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With the PDF defined above, we form the hkehhood function 

£ = l[p{xi), (35) 

i 

where i is an index over events. We perform an unbinned-maximum-hkehhood fit and 
determine the 26 Dahtz plot parameters using the hkehhood function with the signal fraction 
and the lineshape parameters obtained in Sec. 2 and Sec. 3, respectively. 

5. FIT RESULT 

A fit to the 2,824 events in the signal region yields the result listed in Table I. The 
correlation matrix of the 26 parameters after combining statistical and systematic errors 
is shown in the appendix B. Figure 7 shows the projections of the square Dalitz plot in 
data with the fit result superimposed. We also show the mass and helicity distribution for 
each fnr enhanced region along with projections of the fit (Fig. 8). Figure 9 shows the At 
distributions and background subtracted asymmetries. We define the asymmetry in each 
At bin by (iV<j,,g=+i - Nq^,^=_i) / {Nq^^^=+i + Nq^^^=_i), where 7Vg,^g=+i(_i) corresponds the 
background subtracted number of events with q'tag = +1(~1)- The p~7r'^ enhanced region 
shows a significant asymmetry, corresponding to a non-zero value of UZ. 

5-1. Treatment of statistical errors 

With a Toy MC study, we check the pull distribution, where the pull is defined as the 
residual divided by the MINOS error. Here, the MINOS error, which corresponds to the 
deviation from the best fit parameter when —2ln{C/Craa.x) is changed by one, is an estimate 
of the statistical error. Although the pull is expected to follow a Gaussian distribution with 
unit width, we find that the width of the pull distribution tends to be significantly larger 
than one for the interference terms due to small statistics. To restore the pull width to 
unity, we multiply the MINOS errors of the interference terms by a factor of 1.19, which 
is the average pull width for the interference terms obtained above, and quote the results 
as the statistical errors. For the non-interfering terms, we quote the MINOS errors without 
the correction factor. 
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FIG. 7: Distributions of 0' (upper) and m' (lower) with fit results. The histograms are cumulative. 
Solid, dot-dashed, dotted and dashed hatched histograms correspond to correctly reconstructed 
signal, SCF, BB, and continuum PDFs, respectively. 
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FIG. 8: Mass (upper) and helicity (lower) distribution of p^vr (left), p vr^ (middle), and 
(right) enhanced regions. The histograms are cumulative. Solid, dot-dashed, dotted and dashed 
hatched histograms correspond to correctly reconstructed signal, SCF, BB^ and continuum PDFs, 
respectively. 




FIG. 9: Proper time distributions of good tag (r > 0.5) regions for /tag = B^ (upper) and /tag = B^ 
(middle upper), in p+Tr" (left), p~7r~^ (middle), p^n^ (right) enhanced regions, where solid (red), 
dotted, and dashed curves correspond to signal, continuum, and BB PDFs. The middle lower and 
lower plots show the background subtracted asymmetries in the good tag (r > 0.5) and poor tag 
(r < 0.5) regions, respectively. The significant asymmetry in the p~Tr^ enhanced region (middle) 
corresponds to a non-zero value of UZ ■ 
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TABLE I: Results of the time-dependent Dalitz fit. 





Fit Result 


u+ 

TT+, RC 

jT+,Re 
<^-0 

TT+M 


+1 (fixed) 
+1.28 ± 0.13(stat.) ± 0.08(syst.) 
+0.30 ± 0.06(stat.) ± 0.05(syst.) 
+0.62 ± 0.80(stat.) ± 0.57(syst.) 

+0.41 ± 0.52(stat.) ± 0.46(syst.) 
+0.49 ± 0.65(stat.) ± 0.44(syst.) 
+0.86 ± 0.83(stat.) ± 0.49(syst.) 
-0.53 ± 0.39(stat.) ± 0.47(syst.) 
-1.72 ± 0.69(stat.) ± 0.53(syst.) 


UZ 
Uo 

rr-, Rc 
TT-,Re 

TJ-M 


+0.22 ± 0.15(stat.) ± 0.10(syst.) 
-0.62 ± 0.17(stat.) ± 0.09(syst.) 
+0.14 ± 0.11(stat.) ± 0.09(syst.) 
-1.70 ± 1.59(stat.) ± 0.77(syst.) 
-2.46 ± 1.39(stat.) ± 0.86(syst.) 
-0.70 ± 1.59(stat.) ± 0.86(syst.) 
-2.21 ± 1.71(stat.) ± 1.03(syst.) 
-0.83 ± 0.98(stat.) ± 0.65(syst.) 
-0.79+ 1.59(stat.) ± 1.05(syst.) 


1+ 
J_ 

jRe_ 

jRc 
-'+0 
7-Re 
-'-0 
rim 
-'+- 
rim 
''+0 
rim 
-'-0 


-0.03 ± 0.11(stat.) ± 0.06(syst.) 
+0.11 ± 0.11(stat.) ± 0.05(syst.) 
+0.02 ± 0.09(stat.) ± 0.06(syst.) 

+1.62 ± 2.65(stat.) ± 1.23(syst.) 
+1.45 ± 2.41(stat.) ± 1.12(syst.) 
-0.65 ± 1.63(stat.) ± 1.49(syst.) 
-1.76 ± 2.42(stat.) ± 1.31(syst.) 
+0.00 ± 2.06(stat.) ± 1.15(syst.) 
-2.58+ 1.72(stat.) ± 1.33(syst.) 



17 



6. SYSTEMATIC UNCERTAINTIES 



Tables II-IV list the systematic errors for the 26 time-dependent Dalitz plot parameters. 
The total systematic error is obtained by adding each source of systematic uncertainty in 
quadrature. 

6-1. Radial excitations {p' and p") 

The largest contribution for the interference terms tends to come from radial excitations. 
The systematic error related to the radial excitations (p(1450) and p(1700), or p' and p") 
can be categorized into three classes: 1) uncertainties coming from the lineshape variation, 
i.e., the lineshape difference between each decay amplitude, 2) uncertainties in external 
parameters, mp(i45o), rp(i45o), '"^^(itoo)! rp(i7oo), and 3) uncertainties in the common lineshape 
parameters (3 and 7 used for the nominal fit. 

In our nominal fit, we assume all of 6 decay amplitudes have the same contribution from 
p(1450) and p(1700), i.e., we assume equation (15). This assumption, however, is not well 
grounded. In general, the contributions from p(1450) and p(1700) can be different for each 
of the decay amplitudes and thus the systematic uncertainty from this assumption must be 
addressed. Without the assumption about the higher resonances, equation (15) becomes 

= T^fI\s.) , (36) 

where 

F^(s) ^ 5Pr,f,o)(5) + (/5 + AVji^W^pf450)(^) + (7 + AV«)5%f,oo)(5) • (37) 

The variation of the contributions from radial excitations is described by non-zero and 
which are 12 complex variables. We generate various Toy MC samples, where the 
input A'^ and A'^ are fixed but the values of A^/^^ and A*7k arc randomly varied according 
to the constraints on A*/?^ and A''jl; these constraints are obtained from the results in 
Sec. 3, and are combined with the isospin relation [7, 8], which improves the constraints. 
The statistics for each pseudo-experiment are set to be large enough so that the statistical 
uncertainty is neghgible. We assign the variations and the biases of the fit results due to 
the A and A*7k variation as systematic errors. 

For the masses and widths of the p(1450) and p(1700), we quote the values from the 
PDG [19]. To estimate the systematic error coming from uncertainties in their parameters, 
we generate Toy MC varying the input masses and widths and fit them with the masses 
and widths of the nominal fit. Here, we vary the masses by two times the PDG error 
(±50MeV/c^ for p(1450) and ±40MeV/c^ for p(1700)) since the variations between inde- 
pendent experiments are much larger than the la PDG errors, while we vary the widths by 
the ±l(j PDG errors. We quote the mean shift of the Toy MC ensemble as the systematic 
errors. We also take into account the systematic errors from the uncertainties in /3 and 7 
for the nominal fit (equation (30)) in the same way. 
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6-2. SCF 



Systematic errors due to SCF are dominated by the uncertainty in the difference between 
data and MC; these errors are determined from the B ^ D'^*^ p control samples that contain 
a single 7r° in the final state. We vary the amount of SCF by its la error, which is ±100% 
for the CR SCF and +3^°^% the NR SCF in the DS-I(DS-II); we quote the difference 
from the nominal fit as the systematic error. The event fraction for each r region (for CR 
and NR), the wrong tag fractions (for CR) and lifetime used in the At PDF (for CR), which 
are obtained from MC, are also varied and the difference in the fit results is assigned as 
systematic error. 

6-3. Signal Dalitz PDF 

Systematic errors due to the Dalitz PDF for signal is mainly from the Dalitz plot depen- 
dent efficiency. We take account of MC statistics in the efficiency and uncertainty in the 
7r° momentum dependent efficiency correction, e'{p-,^o), obtained from the control samples 
of the decay modes W p-D^*^+, W 7f-D^*^+, ^ p'D^*^'^ and B' tt'D^*^'^. 
The Dalitz plot efficiency obtained from MC is found to have a small charge asymmetry 
(~ 3% at most). We use this asymmetric efficiency for our nominal fit. To estimate the 
systematic error from the asymmetry, we fit the data using a symmetric efficiency and quote 
twice the difference between symmetric and asymmetric efficiencies as the systematic error. 
The Dalitz plot efficiency is r region dependent and obtained as a product with the event 
fraction in the corresponding region, J-'^^^-e^m', 9'), using MC. The difference in the fraction 
for data and MC is estimated to be ~ 10% using the B^ — > D*~tt~^ decay mode as a control 
sample. The fractions are varied by ±10% to estimate the systematic error. 

6-4. Background Dalitz PDF 

The Dalitz plot for continuum background has an uncertainty due to the limited statistics 
of the sideband events, which we use to model the PDF. We estimate the uncertainty by 
performing a Toy MC study of sideband events. With each pseudo-experiment, we model 
the PDF in the same way as we do for real data; with the PDF, we fit the data in the 
signal region; and we quote the variation of fit results as the systematic error. The flavor- 
asymmetry parameters of the continuum background, which are fitted from sideband events, 
are varied by their uncertainties. Systematic uncertainty from the statistics of the BB MC, 
which is used to model the BB Dahtz plot PDF, is also taken into account. 

6-5. TT+TT^TT^ processes other than (P^)° 

The primary contribution to the systematic errors of the non-interfering parameters tends 
to come from the B^ tt+tt^tt^ decay processes that are not B^ (p^)°- We take 
account of the contributions from 5° ^ /o(980)7r°, 5° ^ /o(600)7r°, B° ujn^, and 
non-resonant B^ — > tt'^tt'tt^. Upper limits on their contributions are determined from 
data, except for 5° un^, for which we use the world average of B{B^ — > a;7r°) [20] and 
B{u! TT"*"?!") [19]. For the mass and width parameters of the /o(600) resonance, we use 
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recent measurements by BES [21], CLEO [22], and E791 [23] and take the largest variation. 
We find no significant signals for any of the above decay modes. Using the la upper limits 
as input, we generate Toy MC for each mode with the interference between the — > (p7r)° 
and the other S° — > 7r'*'7r~7r° mode taken into account. We obtain the systematic error 
by fitting the Toy MC assuming B'^ — > (p7r)° only in the PDF. Within physically allowed 
regions, we vary the CP violation parameters of the other 5° 7r"'"7r~7r° modes and the 
relative phase difference between B^ — > (p7r)*^ and the other B^ tt'^tt'tt^ as the input 
parameters, and use the largest deviation as the systematic error for each decay mode. 

6-6. Background fraction 

Systematic errors due to the event-by-event AE-Mi,^ background fractions arc studied 
by varying the PDF shape parameters and the fraction of continuum background, and the 
correction factor to the signal PDF shape by ilex. We also vary the fractions of the BB 
background, which are estimated with MC, by ±50% (±20%) for the 6 — > c (6 — > li) 
transition. 

6-7. Physics parameters 

We use the world average [19, 20] for the following physics parameters: r^o and Anid (used 
for signal and BB background At), the CKM angles of (pi and 02 (used in BB background), 
and the branching fractions oib ^ u decay modes (used in BB background). The systematic 
error is assigned by varying them by ±1(7. The charge asymmetry of B^ — afn^. for which 
there is no measurement and we use zero in the nominal fit, is varied in the physically 
allowed region, i.e., ±1. 

6-8. Background At PDF 

Systematic errors from uncertainties in the background At shapes for both continuum 
and BB backgrounds are estimated by varying each parameter by ±1(7. 

6-9. Vertex reconstruction 

To determine the systematic error that arises from uncertainties in the vertex recon- 
struction, the track and vertex selection criteria are varied to search for possible systematic 
biases. Systematic error due to the IP constraint in the vertex reconstruction is estimated 
by varying the smearing used to account for the B flight length by ±10 iim. 

6-10. Resolution function for the At PDF 

Systematic errors due to uncertainties in the resolution function are estimated by varying 
each resolution parameter obtained from data (MC) by ±la (±2(7). Systematic errors due 
to uncertainties in the wrong tag fractions are also studied by varying the wrong tag fraction 
individually for each r region. 
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6-11. Fit bias 



We observed fit bias due to small statistics for some of the fitted parameters. Since this 
bias is much smaller than the statistical error, we take it into account in the systematic 
errors. We estimate the size of the fit bias by Toy MC study and quote the bias as the 
systematic errors. We also confirm that the bias is consistent between Toy MC and full 
detector MC simulation. 

6-12. Tag-side interference 

Finally, we investigate the effects of tag-side interference (TSI), which is the interference 
between CKM-favored and CKM-suppressed B ^ D transitions in the /tag final state [24]. 
A small correction to the PDF for the signal distribution arises from the interference. We 
estimate the size of the correction using the — > D*~i^i> sample. We then generate MC 
pseudo-experiments and make an ensemble test to obtain the systematic biases. 

TABLE II: Table of systematic errors (1). The notation "< 0.01" means that the value is small 
and less than 0.01, and thus invisible in the number of significant digits shown here. We calculate 
the total systematic error including these small contributions. 
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Other TTTTTT 
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Physics 
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BG At 
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01 
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01 
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Vertexing 
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16 





11 
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Resolution 
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01 


< 





01 
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07 
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03 
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Flavor tagging 


< 


01 
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01 


< 


01 


< 


01 





01 


< 


01 


< 


01 





01 


Fit bias 





01 
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16 





22 





07 





09 





22 





24 


TSI 


< 


01 


< 





01 





01 





01 





01 





03 





01 





01 


Total 


0.08 




0.05 


0.57 


0.46 


0.44 


0.49 


0.47 


0.53 
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TABLE III: Table of systematic errors (2). The notation "< 0.01" means that the value is small 
and less than 0.01, and thus invisible in the number of significant digits shown here. We calculate 
the total systematic error including these small contributions. 
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32 
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Physics 
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02 





06 
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Vertexing 
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03 
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03 
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23 
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05 
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03 





02 





06 





01 


Total 


0.10 


0.09 




0.09 


0.77 


0.86 


0.86 


1.03 


0.65 


1.05 



TABLE IV: Tabic of systematic errors (3). The notation "< 0.01" means that the value is small 
and less than 0.01, and thus invisible in the number of significant digits shown here. We calculate 
the total systematic error including these small contributions. 
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-'-0 


rim 

-'+- 


rim 


rim 
-'-0 


p' and p" 




0.03 




0.02 




0.04 


0.95 


0.59 


1.32 


0.89 


0.84 


0.89 


SCF 




0.01 
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0.01 




0.01 




0.01 


0.33 


0.29 


0.30 


0.31 


0.35 


0.31 


BG Dalitz 




0.01 




0.01 




0.01 


0.34 


0.38 


0.30 


0.32 


0.34 


0.33 


Other TTTTTT 




0.03 




0.03 




0.02 


0.17 


0.15 


0.18 


0.22 


0.15 


0.20 


BG frac. 




0.02 




0.01 




0.01 


0.44 


0.34 


0.33 


0.32 


0.37 


0.29 


Physics 




0.01 




0.01 


< 


0.01 


0.05 


0.06 


0.03 


0.05 


0.05 


0.05 


BG At 


< 


0.01 


< 


0.01 


< 


0.01 


0.05 


0.04 


0.04 


0.05 


0.04 


0.11 


Vertexing 




0.02 




0.02 




0.04 


0.16 


0.28 


0.14 


0.42 


0.37 


0.28 


Resolution 




0.01 




0.01 




0.01 


0.30 


0.21 


0.18 


0.35 


0.25 


0.28 


Flavor tagging 


< 


0.01 


< 


0.01 


< 


0.01 


0.04 


0.07 


0.04 


0.03 


0.07 


0.03 


Fit bias 


< 


0.01 




0.01 


< 


0.01 


0.12 


0.01 


0.27 


0.09 


0.09 


0.22 


TSI 




0.01 


< 


0.01 




0.01 


0.09 


0.07 


0.04 


0.04 


0.03 


0.07 


Total 




0.06 




0.05 




0.06 


1.23 


1.12 


1.49 


1.31 


1.15 


1.33 
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TABLE V: Correlation matrix of the quasi-two-body parameters, with statistical and systematic 
error combined. 







c 


AC 


s 


A5 




+1.00 










c 


-0.15 


+1.00 








AC 


+0.08 


+0.19 


+1.00 






S 


+0.02 


-0.02 


-0.00 


+1.00 




AS 


-0.02 


-0.00 


-0.01 


+0.30 


+1.00 



A^^ - 


-0 


12 + 0.05 + 0.03 , 


c = 


-0 


13 + 0.09 + 0.06 , 


AC = 


+0 


35 + 0.10 + 0.06 , 


S = 


+0 


06 ±0.13 ±0.07, 


AS = 


-0 


12 ±0.14 ±0.07, 



7. QUASI-TWO-BODY PARAMETERS 

We calculate quasi-two-body CP violation parameters as 

C.^, AC.^, A..^. (39) 

We obtain 

(40) 
(41) 
(42) 
(43) 
(44) 

where first errors and second errors are statistical and systematic, respectively. The corre- 
lation matrix is shown in Table V. 

One can transform the parameters into the direct CP violation parameters Ap~ and A~J~ 
defined as 

ACP , /J , ACPAr 

'^p- - 1 + AC + ^^/C ' ^ > 

A^^ -C- A^^AC 

A-+ = 1^ I ^ f46l 

^p- - i-AC-A^^C ' ^ ^ 

which can be more convenient for interpretation. We obtain 

A+- = +0.22 ± 0.08 ± 0.05 , (47) 
= +0.08 ±0.17 ±0.12 , (48) 

with a correlation coefficient of +0.53. Our result differs from the case with no direct CP 
asymmetry {Ap~ = and A'^ = 0) by 2.4 standard deviations (Fig. 10). Measurements 
with a larger data sample may be needed to observe direct CP violation. 
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FIG. 10: Contour plot of the confidence level for the direct CP violation parameters vs. AfJ~. 



The ratio of the branching fractions B{B^ p^n^'*)/B{B^ — > p^n^) can be related to our 
measurement as 



Our measurement yields 



B{B^ p±7rT) 



0.133 ±0.022 ±0.023 . 



(49) 



(50) 



This is consistent with the dedicated quasi-two-body branching fraction measurement of 
50 ^ ^0^0 f^Q^ Bg^g [25]: 



0^0\Belle 3 -|^2+0-88+0.60 



0.82-0.76 _ n. -, qn+0-049 
— — — — — U.10U_o.o46 > 



23(50 ^ p±7rT)WA 24.0 ± 2.5 

where we use the world average [20] for the denominator. 

We also measure the CP violating parameters of B^ —>■ p°7r° calculated as 

A ^0 . 2/0 



We obtain 



(51) 



(52) 



Apo^o = -0.45 ± 0.35 ± 0.32 , 
Spo^o = +0.15 ± 0.57 ± 0.43 , 



(53) 
(54) 



with correlation coefficient of —0.07. Our measurement of ApO^^o is consistent with the 
previous measurement from Belle [25] . This is the first measurement of SpO^^o . 
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8. CONSTRAINT ON THE CKM ANGLE 02 



We constrain the CKM angle 02 from our analysis following the procedure described in 
Ref. [6]. With three 5° {p^T decay modes, we have 9 free parameters including 

9 = (6 complex amplitudes = 12 d.o.f.) + 02 /rr\ 
— (1 global phase) — (1 global normalization) — (1 isospin relation = 2 d.o.f.) , 

where we make use of an isospin relation that relates neutral B decay processes only [7, 8]. 
Parameterizing the 6 complex amplitudes with 9 free parameters, we form a function 
using the 26 measurements of our time-dependent Dalitz plot analysis as constraints. We 
first optimize all the 9 parameters to obtain a minimum x^, Xmin! "^^ then scan 02 from 0° to 
180° optimizing the other 8 parameters, whose resultant minima are defined as X^(02); and 
Ax^(02) is defined as Ax^(02) = X^(02) ~ Xmm- Performing a Toy MC study following the 
procedure described in Ref. [26], we obtain the 1 — C.L. plot in Fig. 11 from the Ax^(02)[28]. 

In addition to the 26 observables obtained from our time-dependent Dalitz plot analysis, 
we use the following world average branching fractions and asymmetries: B{B° — > p^7r^), 
B{B+ p+ttO), -> p+ttO), B{B+ p07r+), and A{B+ p^T^+) [20], which are not 

correlated with our 26 observables. With the 31 measurements above, we perform a full 
combined Dalitz and isospin(pentagon) analysis. Having 5 related decay modes, we have 12 
free parameters including 02: 

12 = (10 complex amplitudes = 20 d.o.f.) + 02 

— (1 global phase) — (4 isospin relations = 8 d.o.f.) . 

The detail of formation can be found in appendix C. The x^^j^ obtained is 9.8, which is 
reasonable for 31 (measurements) — 12(free parameters) = 19 degrees of freedom. Following 
the same procedure as above, we obtain the 1 — C.L. plot in the Fig. 12. We obtain 02 = 
(83^23)° as the central value and la errors (corresponding to 68.3% C.L.). A large CKM- 
disfavored region (02 < 8° and 129° < 02) also remains. 



9. CONCLUSION 

We have performed a full Dalitz plot analysis of the tt+tt^tt'^ decay mode. The 

results are consistent with no direct CP violation and the previous measurement of B^ — > 
p°7r° branching fraction and flavor asymmetry. Combining our analysis and information 
from charged B decay modes, a full Dalitz plot and isospin analysis is performed to obtain 
a constraint on 02. We obtain 02 = (83^23)° as the central value with la errors. However, a 
large CKM-disfavored region (02 < 8° and 129° < 02) also remains. In principle, with more 
data we may be able to remove the additional 02 solutions. 
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FIG. 11: 1 — C.L. vs. 02 obtained by Dalitz analysis. The dotted horizontal line corresponds to 
the la confidence level. 
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FIG. 12: 1 — C.L. vs. (/>2 obtained by a full combined Dalitz and isospin(pentagon) analysis. The 
dotted horizontal line corresponds to the la confidence level. 
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APPENDIX A: PDFS FOR TIME-DEPENDENT DALITZ PLOT ANALYSIS 



In this section, we describe the details of the PDF for each component, which appear in 
equations (31)-(34). 



1. Signal PDF 



a. PDF for correctly reconstructed events 



In terms of the event fractions for the l^^ fiavor tagging region (^true)' Dalitz plot 
dependent efficiency (e'), the tt*^ momentum dependent efficiency correction taking account 

of the difference between data and MC (e'), wrong tag fraction {wi), and the wrong tag 
fraction difference between and (Aw/), the PDF for correctly reconstructed events is 
given by 

Pivnei^E, Mbe; m' , e'; At, gtagiP^ro) = 

P ■ ptrueiAE, Mbc;p,o) • e'(m', 9') e'(p,o) • |det J(m', 9')\ ■ pl^,{m', Ai, feg), ^ ' 

where 



pl,,{m',9';At,q,^,) = ■ {(1 - gtagAt.0(l^3.P + l^s.P) 

+gtag(l - 2wi) ■ [-(I^sttI' - I^ttH cos(AmdAi) + 2Im[M3^Ay sin(AmdAi) 



(A2) 



For the At PDF, the above equation is convolved with the resolution function [14]. The 
A£'-Mbc PDF is normahzed such that 

[[ fiA^c?MbePtrue(A£;,Mbc;p,o) = l (Vp,o) , (A3) 

J J signal region 

since we define the Dalitz plot efficiency for events inside the signal region. With the PDFs 
in the At-gtag direction being also normalized to be unity, the integral inside the signal 
region, ntrue, is, 

nuue = ^true > (A4) 
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^Le = Y. 1^^^ IfdAEdM^, 11 dm' d9' pl^^{AE,M^,;m',9';At,qt,^;p^o) (A5) 

9tag signal region J J SDP, Veto 

= -^true / / dm' d9' ei,,,(m', 9') e'(p,o) |det J| (| + \As^\^) , (A6) 

^ ^SDP, Veto 

where the correlation between p^o and m' is properly taken into account in the integration 
of the last line. The notation f f^^^ ^^^^ dm' d9' means integration over the square Dalitz 
plot with the vetoed region in the Dalitz plot taken into account. 

The 7r° momentum dependent AE-M\^c PDF, ptrue(A£^, Mbc;p7ro), is modeled using MC- 
simulated events in a binned histogram interpolated in the p-^o direction, to which a small 
correction obtained with B° — > p~ is applied to account for the difference between MC 
and data. 

The Dalitz plot distribution is smeared and distorted by detection efficiencies and de- 
tector resolutions. We obtain the signal Dalitz plot efficiency from MC to take the former 
into account. We introduce a dependence of the efficiency on the r region, e^jue) since a 
significanct dependence is observed in MC. Small corrections, e'(p^o), arc also applied to 
the MC-determined efficiency to account for differences between MC and data. We use 
W p-D^*^+, W TT'D*-, B- p-D^ and B' t^'D^ decays to obtain the correc- 
tion factors. The detector resolutions are small compared to the widths of p(770) resonances; 
this is confirmed by MC to be a negligibly small effect. 



h. PDF for SCF events 

Approximately 20% of signal candidates are SCFs, which are subdivided into ~ 4% of 
NR SCF and ~ 16% of CR SCF. It is therefore important to develop a model that describes 
the SCF well. The time-dependent PDF for SCF events is defined as 

p\{AE, Mbe; m', 9'- At, g^ag) = -Pii^E, My,,; s^) -p^im', 9'; At, gtag) , (^ = NR, CR) (A7) 

where is the event fraction in each tagging r-bin and 

p,{m', 9'; At, g^ag) = ■ {(1 - gtagA^,^) p\^\n^ , 9') 

+gtag(l - 2w,^) • [-pf°"(m', 9') cos(AmdAi) +pf°(m', 9') sin(AmdAi)]| . 

(A8) 

The AE-M\,c PDF is normalized inside the signal region as 



1 1 dAE dMy,, pi{AE, Mbc; Si) = 1 (Vs^) . (A9) 

J J signal region 

With the PDFs in the At-q^ag direction being normahzed to unity and Ylii^l — 1> 
integral inside the signal region, nj, is 

= J] 5^ /(iAt / jdAE dMy,, j I dm' d9' p\{AE, Mbc; m' , 9'; At, gtag) (AlO) 

^ J J J signal region J J SDP, Veto 



gtag 



dm' d9' pY'%m', 9') . (All) 

SDP, Veto 
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We find that the A£'-Mbc distribution for SCF has a sizable correlation with Dalitz 
plot variables, but with only one of its two dimensions. We thus introduce a model with 
dependences on the Dalitz plot variable Sj. The variable scr = s± = max(s+,s_) is used, 
because the CR SCF can be divided into a 7r+ replaced SCF and a tt" replaced SCF, where 
s_ is used for tt"*" (tt") replaced SCF. Here, wc exploit the fact that almost all of 
the TT"*" (tt") replaced SCF distributes in the region of s+ < s_ (5+ > s_). For the NR 
SCF, snr = Sq. This parameterization models the correlation quite well, with each of the 
parameters Si reasonably related to the kinematics of replaced tracks. 

Since track (tt) replacement changes the measured kinematic variables, the SCF events 
"migrate" in the Dalitz plot from the correct (or generated) position to the observed posi- 
tion. Using MC, we determine resolution functions -Rj(^obs> ^obs5 "^gcni ^gcn) describe this 
"migration" effect, where {m'^^^, 9'^^^) and {m'^^^, 9'^^^) are the observed and the generated 
(correct) positions in the Dahtz plot, respectively. Together with the efficiency function 
^i("^gen> ^gen)> which is also obtained with MC, the Dalitz plot PDF for SCF is described as 

p>{m', 9') = [{R. ■ e.) ® p^,,] (m', 9') (A12) 

dm'^en d^'gen , ^' ^gen' ^gen) " ^^("^gen' ^gen) " Pgen("^gen' ^gen) (^1^) 

SDP 



where 



( i = CR, NR , J = Life, Cos, Sin ) 



Pgen«en, 0',J = \detJ\{\As,\' + \As^\') , (Al4) 

P?e°„^Kcn, C) = |det J|(|A3.|' - |:43.n , (Al5) 



Pi^Hen,Ogen) = |detJ|2Im 



- 4o A* 



(A16) 



For the NR SCF, the shape of the At PDF defined in equation (A8) is exactly the same 

as correctly reconstructed signal, i.e., tnr = r^o, wi^^ = wi, and Awi^^ = Awi, since the 
replaced track, 7i^\ is not used for cither vertexing or flavor tagging. On the other hand, for 
the CR SCF, the At PDF is different from correctly reconstructed signal, since the replaced 
is used for both vertexing and flavor tagging. Thus, we use MC-simulated CR SCF events 
to obtain tqr, wp^, and Awi'"^, which are different from those of correctly reconstructed 
signal events. In particular, Awi'"^ is opposite in sign for the Tr~^ and vr" replaced SCFs, 
which is due to the fact that the replaced tt^ tends to be directly used for flavor tagging in 
the slow pion category. 

2. Continuum PDF 

The PDF for the continuum background is 
p^gg{AE,My,^;m',9';At,qt^s) 

= • p[,(AE, Mb.) • p,,{m', 9'- AE, M^J • [i±SiM|K^j . p^^^At) , ^^^^> 

where J^g^ is the event fraction for each r region obtained in the signal yield fit. All the 
terms on the right hand side of the equation are normalized to be unity so that 

I ^^bc / / dm' d9' p[g{AE, Mbe; m', 9'; At, q,,^) = 1 . (A18) 

I 9tag signal region J J SDP, Veto 
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Since the allowed kinematic region is dependent on AE" and Afbc, the Dalitz plot distri- 
bution is dependent on AE and Mbc- We define a AE-My^c independent PDF, Pqq{rn'^ca.iey 
where m'^^aie ^ re-defined SDP variable with the kinematic effect taken into account as 

m' ,^ = - arccos 2 -. ^ 1 , A19) 

where 

AMbc = Mbc - . (A20) 
Using the A£^-Mbc independent PDF, pqq{m', 9'; AE, Mbc) is described as 

p„(m', 0; AE. M,J = ^(^^^^^^^ ■ ■ 0') (A21) 

for the region, m™''^ < mo < min(m™'^^, m™*^^ + AE + AMbc) {Pqq = otherwise), where 
N{AE + AMbc) and sin(7rm')/ sin(7rmg^^jp) are a normalization factor and the Jacobian for 
the parameter transformation of m'^caie ~^ respectively. We obtain the Pqqi'fn'^ca.ie^ 
distribution from data in part of the sideband region, — O.lGeV < AE < 0.2 GeV and 
5.2 GeV/c^ < Mbc < 5.26 GeV/c^ where the contribut ion from BB background is negligible. 

Since we find significant fiavor asymmetry depending on the location in the Dalitz plot, 
we introduce the following term to take account of it: 

which is r region dependent. The asymmetry is anti-symmetric in the direction of 9' , i.e., 
A\m\9') > {A^{m',9') < 0) in the region of 9' > 0.5 {9 < 0.5), and the size of the 
asymmetry is ~ 20% at most in the best r region. This effect is due to the jet-like topology 
of continuum events; when an event has a high momentum vr" (tt"*") on the CP side, the 
highest momentum vr on the tag side tends to have + (— ) charge. The highest momentum 
TT on the tag side with -|- (— ) charge tags the fiavor as 5° (B^)- Since an event with a high 
momentum vr" (7r+) resides in the region 9' > 0.5 {9' < 0.5), a continuum event in the region 
9' > 0.5 {9' < 0.5) tends to be tagged as S° (B^). We again parameterize the A\m',9') in 
a AE'-Mbc independent way as 

A\m',9')^A\m',,,,,,9'), (A23) 

and model with a two-dimensional polynomial, whose coefficients are determined by a fit to 
data in the sideband region. 



3. BB background PDF 

The treatment of BB background is different for CP eigenstate modes and fiavor specific 
or charged modes. The PDF for the CP eigenstate modes is 

p'bb{^E, Mbc; m', 9'; At, feg) = ^bb • Pbb{AE, Mbc) • PBB{m', 9') ■ pBB{At, gtag) , (A24) 

where pBB{At, ^tag) is a time-dependent CP violation PDF normalized as 

^ j dAt pBsiAt, gtag) = 1 • (A25) 

9tag 
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For the flavor specific or charged modes, the PDF is 

p^Ssi^E, Mbc; m', 9'; At, gtag) = ^BB -pBBiAE, Mbc) ^PBB{m', 9'; Qrec) -pBBi^t, gtag, qrec) , 

Qrec 

(A26) 

where the Dahtz plot PDF Pb_b(^'i 9'; grcc) is dependent on the true flavor of the CP (fully 
reconstructed) side, greo and the time dependent part is a mixing PDF (lifetime PDF with 
flavor asymmetry) for flavor specific (charged) modes. The At PDF is normalized as 

/^^* PBB{At, feg, grec) = 1 ■ (A27) 

Q'tag Qrec 

The AE-M^c PDF and Dalitz plot PDF are obtained mode-by-mode from MC. The 
Dahtz plot PDF of the CP eigenstate modes is assumed to have following symmetry 

PBB{m', 9') = PBB{m\ 1-9'), (A28) 

while that of flavor specific and charged modes is assumed to have following symmetry 

PBBim', 9'; q,ec) = PBB{m', 1 - 9'; -g^ec) • (A29) 

The total PDF of BB background is a linear combination of each mode with efficiencies and 
branching fractions taken into account. 

APPENDIX B: CORRELATION MATRIX OF THE FIT RESULT 

Tables VI- VIII show the correlation matrix for the 26 parameters determined in the time- 
dependent Dalitz plot analysis, corresponding to the total error matrix with statistical and 
systematic error matrices combined. We assume no correlation for the systematic errors. 

APPENDIX C: METHOD OF 02 CONSTRAINT 
1. Formalism 

We define amphtudes as 

A+ 
A- 



= A{B' p+TT-) , (CI) 

= AiB'' p-7r+) , (C2) 

= A(S° ^ p°7r°) , (C3) 

^ A{B+ p+7r°) , (C4) 

^ AiB-' p\+) , (C5) 
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TABLE VI: Correlation matrix (1) of the 26 fitted parameters, with statistical and systematic 
errors coml)iiicd. 







K 




^ ^1 Up 


^t4- l-ip 

^-0 


7-7-+,lm 


T7-+,im 


7-7-+,lm 




+1.00 


















+0.21 


+1.00 














uXr 


+0.05 


+0.03 


+1.00 














+0.09 


+0.01 


+0.01 


+ 1.00 










1 


-0.02 


-0.09 


+0.01 


+0.01 


+ 1.00 










+0.02 


+0.01 


+0.03 


+0.00 


-0.00 


+1.00 






rr+,Im 


-0.04 


-0.09 


-0.00 


+0.14 


+0.02 


-0.00 


+1.00 






-0.13 


-0.10 


-0.01 


-0.02 


+0.02 


-0.00 


+0.02 


+ 1.00 




+0.05 


+0.02 


-0.00 


-0.01 


+0.00 


-0.02 


-0.01 


-0.01 


uz 


-0.22 


-0.07 


-0.03 


-0.03 


-0.03 


-0.01 


+0.01 


+0.02 




+0.05 


+0.09 


+0.01 


+0.00 


-0.04 


+0.00 


-0.05 


-0.06 




-0.04 


-0.02 


-0.06 


-0.01 


-0.00 


+0.02 


+0.00 


+0.01 




-0.05 


-0.02 


-0.00 


-0.13 


-0.00 


+0.00 


-0.04 


+0.02 


^ ^ "Dp 


-0.02 


-0.05 


-0.00 


-0.00 


+0.08 


-0.00 


+0.02 


+0.07 




-0.04 


-0.02 


-0.00 


-0.01 


-0.00 


+0.05 


+0.00 


+0.01 


TT-,Im 


-0.04 


-0.08 


-0.01 


-0.06 


+0.01 


-0.00 


-0.08 


+0.02 


[/-^"^ 


-0.02 


-0.02 


-0.00 


-0.00 


-0.03 


-0.00 


+0.00 


-0.12 


^+ 


+0.00 


+0.00 


-0.01 


-0.01 


-0.00 


-0.02 


-0.02 


-0.00 


/_ 


+0.07 


+0.03 


-0.01 


+0.01 


-0.01 


+0.03 


-0.00 


+0.01 


^0 


+0.01 


+0.01 


+0.00 


+0.01 


+0.01 


+0.00 


-0.02 


-0.03 


jRe_ 


-0.02 


-0.00 


+0.01 


+0.00 


+0.00 


-0.20 


+0.00 


+0.00 


-'+0 


+0.00 


+0.02 


+0.00 


-0.12 


-0.01 


+0.00 


-0.01 


+0.00 


rRe 
-'-0 


-0.04 


+0.02 


-0.01 


-0.01 


-0.11 


+0.00 


-0.01 


-0.25 


rim 
-'+- 


-0.02 


-0.01 


+0.09 


-0.00 


-0.00 


+0.03 


+0.00 


+0.00 


rim 
-'+0 


+0.01 


-0.02 


+0.00 


+0.04 


+0.00 


+0.00 


+0.07 


+0.00 


rim 

-'-0 


-0.07 


-0.04 


-0.01 


-0.01 


-0.09 


-0.01 


+0.01 


+0.08 
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TABLE VII: Correlation matrix (2) of the 26 fitted parameters, with statistieal and systematic 
errors coinfnncd. 
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TT—,t\.e 




TT—,n.e 
^-0 










+1.00 


















uz 


-0.06 


+1.00 


















+0.00 


-0.00 


+1.00 
















-0.04 


+0.02 


-0.00 


+1.00 












TT ,«.e 


-0.17 


+0.03 


-0.08 


+0.01 


+1.00 










TT .itC 


+0.01 


-0.10 


-0.18 


-0.00 


+0.02 


+ 1.00 










+0.06 


-0.00 


-0.00 


+0.16 


-0.01 


+0.00 


+1.00 








-0.00 


+0.01 


-0.08 


+0.00 


+0.03 


+0.02 


+0.00 


+1.00 






-0.00 


+0.02 


+0.01 


-0.00 


-0.00 


+0.02 


+0.00 


+0.00 


+ 1.00 


1+ 


-0.02 


-0.00 


-0.00 


+0.02 


+0.00 


+0.00 


-0.02 


+0.01 


+0.01 


/_ 


-0.00 


-0.02 


+0.01 


+0.02 


-0.00 


-0.04 


-0.01 


-0.01 


-0.08 




-0.00 


-0.01 


+0.07 


+0.00 


+0.01 


+0.02 


-0.00 


-0.04 


-0.04 


jRe 


+0.02 


+0.01 


-0.00 


-0.04 


-0.00 


-0.00 


-0.19 


-0.00 


+0.00 


7-Re 
-'+0 


-0.00 


+0.00 


+0.01 


-0.00 


+0.06 


-0.01 


+0.00 


+0.11 


+0.01 


rRc 
-'-0 


-0.01 


+0.07 


+0.05 


+0.01 


-0.00 


-0.09 


+0.00 


-0.01 


+0.13 


rim 


+0.01 


+0.03 


-0.00 


+0.06 


-0.00 


-0.00 


+0.05 


+0.00 


+0.00 


rim 

''+0 


+0.02 


-0.00 


-0.01 


-0.00 


-0.25 


-0.00 


+0.00 


-0.02 


+0.01 


rim 
^-0 


-0.00 


+0.03 


-0.01 


+0.00 


+0.00 


+0.15 


+0.01 


+0.02 


+0.10 



TABLE VIII: Correlation matrix (3) of the 26 fitted parameters, with statistical and systematic 
errors combined. 





1+ 


/_ 




jRe_ 


rRe 
-'+0 


rRe 
-'-0 


rim 

-'+- 


rim 


rim 

-'-0 


1+ 


+1.00 


















7_ 


-0.06 


+1.00 
















lo 


+0.01 


+0.02 


+1.00 














jRe 


-0.08 


-0.00 


-0.00 


+1.00 












Ifo 


-0.02 


-0.00 


-0.10 


+0.00 


+1.00 










rRe 


-0.01 


+0.14 


+0.00 


+0.00 


+0.00 


+ 1.00 








rim 


-0.03 


-0.02 


-0.00 


-0.28 


+0.00 


+0.00 


+1.00 






rim 

''+0 


-0.13 


+0.01 


-0.12 


+0.01 


-0.11 


-0.00 


+0.01 


+1.00 




rim 
^-0 


+0.01 


-0.12 


-0.19 


+0.00 


+0.03 


+0.03 


+0.00 


+0.03 


+ 1.00 
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and 

(C6) 
(C7) 
(C8) 
(C9) 

These amplitudes are obtained from 1) 26 measurements determined in the time-dependent 
Dahtz plot analysis as well as 2) branching fractions and asymmetry measurements, and 
give a constraint on 02- 

Equations (21)-(25) define the relations between the amplitudes for the neutral modes 
and the parameters determined in the time-dependent Dalitz plot analysis. The relations 
between the branching fractions and asymmetries, and the amplitudes are 

B(p±7r^) = c ■ + + + |:4-|) • tbo , (Cll) 

B(p+7r°) = c . + \A-r) ■ , (C12) 

B(p°7r+) = c • + \A'-\') ■ Ts+ , (C13) 

14-012 _ I 4+012 

Ap'-^) = '|^o-|2 + |^o+j2 ' (CIS) 

where c is a constant and the lifetimes r^o and Tb+ are introduced to take account of the total 
width difference between B^ and B^. Note that we do not use quasi- two-body parameters 
related to neutral modes except for S(p^7r^), since they are included in the Dalitz plot 
parameters. 

The amplitudes are expected to follow SU (2) isospin symmetry to a good approximation 

A+ + A- + 2A'' = A+ + A- + 2A^ ^ 

= ^(A+o + A^+) = V2{A-^ + A"-) , ^ ^ 

^+0 _ ^0+ _ y2(A+ - A-) = I-° - - V2{A- - 1+) , (C17) 

where _ _ _ 

A" = e-^^-^M^ ^ ^-0 = e-2i<^2^-o ^ ^0- = e-2i02^o- _ ^^^g^. 

Note that there is an inconsistency in equation (CI 7) between Ref. [7] and Ref. [8]; we follow 
the treatment of Ref. [7], which we believe is correct. 



2. Parameterization 

Here we give two examples of the parameterization of the amplitudes. The first example 
may be more intuitive, while the second example is well behaved in the fit. The results are 
independent of the parameterizations with respect to the constraint on 02- 
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a. Amplitude parameterization 



We can parameterize the amplitudes as follows [7] 

A+ = e-^'^2r+ + P+ , (C19) 

A- = e-^'I'^T- + P- , (C20) 

^0 ^ g-i02yo _ l(p+ + p-) ^ (C21) 

v/2A+° = e-^'^2T+° + P+ - P- , (C22) 

72^1°+ = e-''^\T+ + T- + 2T° - T+°) - P+ + P" , (C23) 



and 

= e+^-^^T- + P- , (C24) 
A- = e+''^^T+ + P+ , (C25) 

Z° = e+^'^2T° - ^(P+ + P-) , (C26) 

v/2A-° = e+"^''T+° + P+ -P- , (C27) 
V2A°- = e+^^2(r+ + T- + 2r° - T+°) - P+ + P" , (C28) 

where the overall phase is fixed with the convention ImT+ = 0. Thus, there are 6 complex 
amplitudes, T'^ ,T' , P^ , P~ , and T+°, corresponding to 11 degrees of freedom; and 02, 

corresponding to 12 degrees of freedom in total. This parameterization automatically sat- 
isfies the isospin relations without losing generality, i.e., the isospin relations are the only 
assumption here. 



b. Geometric parameterization 



We can parameterize the amplitudes using the geometric arrangement of Fig. 13 that 
satisfies the isospin relation of equation (C16). This figure is equivalent to Fig. 3 of Ref. [8], 
except that the sides corresponding to P° — > p~Ti^ and P° p^-K^ are swapped. This 
difference is not physically significant. We apply this modification only to obtain a better 
behaved parameterization; the parameterization here uses the angles a;_ and 6_ related to 
the process P°(P°) — > p'n'^, which are better behaved than those related to P°(P°) — > p°7r°. 

To parameterize the amplitudes, we use 4>2 and the following 11 geometric parameters: 

cu', 9^, 6+, 6_, b', a+, a_, L, (C29) 
where b and a imply branching fraction and asymmetry, respectively. In terms of these 



35 



V2A' 

FIG. 13: Complex pentagons formed from the B ^ pir decay amplitudes. 



parameters, the amplitudes can be described as follows 



A+ = e'('^++^+/2V&+(l - a+)/2 , (C30) 

X+ ^ e'('^+-^+/2V^'+(l + «+)/2 , (C31) 

A- = e'('^-+^-/'V^'-(l-«-)/2, (C32) 

I- = e'('^--^-/'V^'-(l + «-)/2 , (C33) 

= (L->1+ - A-)/2 , (C34) 

1° = (L-I+ -I-)/2 , (C35) 

= e'^^'^/WJ^, (C36) 

= ^ - A«+ , (C37) 

■^+0 _ ^0+ _ ^(^+ _ ^ _ /2 , (C38) 



2\/2 

= ^ - I°- . (C39) 

Equation (C38) exploits the isospin relation of equation (C17), which Fig. 13 does not incor- 
porate geometrically. The phase 02 enters when the As are converted into A's with equation 
(C18). When we perform the analysis only with the time-dependent Dahtz plot observables 
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and without the information from charged decay modes, we remove the parameters to' and 
6' from the fit and fix L to be a constant. 

This geometric parameterization has a substantial advantage in terms of required com- 
putational resources, compared to the parameterization based on the T and P amplitudes 
described in the previous section. In the procedure to constrain 02) the minimum ^ has to 
be calculated for each value of 02- To avoid local minima, initial values of the parameters for 
the minimization have to be scanned and this inflates the computing time, which increases 
exponentially with the number of parameters. However, the number of parameters to be 
scanned decreases in the geometric parameterization. Among 11 parameters except for 02, 
five of them, 6+, 6-, 6', a+, and a_, are related to the branching fractions and asymmetries. 
Since in most cases they do not have multiple solutions, we do not have to scan the initial 
values of them. In addition, the optimum initial value for L can also be determined using 
other parameters and 6o, the nominal branching fraction of — > p°7r°, from the following 
relation 

6o = I - e^"+ - e^"- aA72 f , (C40) 

up to a two-fold ambiguity. Here feo is calculated using the input parameters as 

^0^^^-^^, (C41) 

based on equations (21) and (Cll). The explicit solution for the optimal initial value of L 
is 

L = Re7±-^6o-(Im7)' ( = e'""^ ^ e'''- ^hlj2 ) . (C42) 

When 6o ~ (Ij^t)^ < 0, there is no real- valued solution and L — Re7 is the optimum initial 
value. With the optimum values calculated above, the initial value of L does not have to 
be scanned, except for the two fold ambiguity. Consequently, the number of parameters 
to be scanned in this parameterization is only five, corresponding to a;+, a;_, a;', 6*+, and 
while all of 11 or maybe 10 parameters have to be scanned in the T and P amplitude 
parameterization. This leads to a substantial reduction of the computational resources 
required. 
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